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Abstract 

Poisson brackets for the Hamiltonian dynamics of vortices are dis- 
cussed for 3 regimes, in which the dissipation can be neglected and 
the vortex dynamics is reversible: (i) The superclean regime when the 
spectral flow is suppressed, (ii) The regime when the fermions are 
pinned by crystal lattice. This includes also the regime of the extreme 
spectral flow of fermions in the vortex core: these fermions are ef- 
fectively pinned by the normal component, (iii) The case when the 
vortices are strongly pinned by the normal component. All these lim- 
its are described by the single parameter Co, which physical meaning 
is discussed for superconductors containing several bands of electrons 
and holes. The effect of the Fermi-surface topology on the vortex 
dynamics is also discussed. 

1. Introduction. The problem of the vortex motion is well understood for 
the translational invariant Fermi superfluids Jl], ^ [| |J. When the theory 
is applied to superconductors, one should take into account the effect of 
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the band structure of the crystal on the vortex dynamics. In this case the 
topology of the Fermi- surface is to be important and we make an attempt to 
consider this effect in the regimes, when the dissipation is small and can be 
neglected. 

2. Poisson Brackets Formalism for the Vortex Motion. We start with the 
phenomenological hydrodynamic equations for the system of distributed vor- 
tices in superffuids or superconductors. In the limit of vanishing dissipation 
the dynamics of the collective variables becomes conservative and in principle 
can be described by the effective action. However, as usually occurs in the 
hydrodynamic systems, such action is not well defined and the Hamiltonian 
formalism in terms of the Poisson brackets (PB) becomes preferrable. In 
this formalism the Hamiltonian is the function of the relevant hydrodynamic 
variables characterized by the algebra of the Poisson brackets. 

We assume that the normal component is clamped, ie its velocity v n = 0. 
This is typical for superfluid 3 He due to its high viscosity and for super- 
conductors where v ra is fixed by the impurities in the crystal lattice. The 
remaining hydrodynamic variables at low temperature are the mass density 
p and the superfluid velocity v s , which is non-potential in the presence of the 
distributed vorticity. The Hamiltonian 

H = Jd 3 r[^vy s + e(p)) , (2.1) 

contains the internal energy density e(p) and the kinetic energy. In crystals, 
due to absence of the Galilean invariance, the superfluid component does not 
coincide with the density p of the electrons even at T = 0: p s (T = 0) ^ p. 
Actually p s can be much smaller than p because most of the electrons are 
concentrated in the completely filled bands. Nevertheless the hydrodynamic 
equations are valid for the descrition of the long-wave-length dynamics of p 
and v s . 

The motion equations are obtained as the Liouville equations 

! = {*,,} , £ = Wv.} . (2.2) 

if one uses the PB between the variables. These PB are universal, ie they do 
not depend on the Hamiltonian ||, and we propose: 

{p(r),v s (r')} = W(r-r') , (2.3) 
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{p(r),p(r')} = , (2.4) 

K(r), ^-(r')} = -e ljk {WXV n s)k 6(r - r') . (2.5) 

P — 

The first two PB are conventional (see jS|), eg the Eq.(2.3) follows from the 
fact that the particle number and the phase of the condensate are canonically 
conjugated variables. The Eq.(2.5) contains a new variable Co, which is to 
be the dynamical invariant of the system, = 0. With this constraint the 
Poisson brackets satisfy the Jacobi identity {a{bc}} + {b{ca}} + {c{ab}} = 0. 
The Eq.(2.5) with Co = was derived in superfluids in the T = limit |J. 
In derivation it was assumed that each element of the vortex moves with the 
local superfluid velocity v s (r), which corresponds to the Helmholtz theorem 
for the perfect liquid. The PB Eq.(2.5) with Co = oo was written in ||. This 
corresponds to the motion of vortices with the local normal velocity v n (r). 
This is actually a true hydrodynamic regime, which requires the complete 
pinning of the vortex lines by the normal component or by the heat bath of 
the crystal lattice JF). 

The range < Co < oo corresponds to the intermediate regimes when the 
vortex is unpinned, but some groups of particles (electrons) are pinned by the 
heat bath during the vortex motion. In addition to the trivial localization 
of the electrons by the crystal lattice and impurities, the quasiparticles are 
also pinned by the anomalous process of the spectral flow in the core of the 
vortex, discussed in || f|. When the vortex moves in the regime of 
the extreme spectral flow, the momentum is effectively transferred from the 
vortex to the heat bath, which corresponds to the pinning of fermions by the 
heat bath. 

With the PB in Eqs. (2.3-5), the Liouville equations (2.2) of the nondissi- 
pative dynamics become 

|| + V-(/3 s v s ) = , (2.6) 

o^v 1 S II 

-£ + ^(p> s )x(Vxv s )+V/i = ,//=—. (2.7) 

at p — C dp 

Using the kinematic definition of local velocity of vortex lines[[[I], |J: 
<9 t v s + Vp = -v L x (V x v.) , (2.8) 
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one obtains the following relation between and v s : 



(p - C )v L = p s v s 



(2.9) 



For superconductors the lhs of Eq.(2.9) determines the Hall conductivity in 
the limit of vanishing dissipation: 



where B is magnetic field. 

The equation (2.9) is not Galilean invariant. The Galilean invariance 
is restored by introducing the velocity of the normal component v n , which 
coincides with the velocity of crystal lattice in the case of superconductors. 
Then multiplying by the circulation k (with \k\ = Nnh/m, and N being the 
winding number of the vortex) one obtains the equation for the balance of 
Magnus, spectral flow and Iordanskii forces acting on the vortex |3J: 



where p n = p — p s is the tensor of the normal density. 

The density p and the superfluid density tensor p s are well determined 
quantity, the latter being determined by the current- current correlation func- 
tion. Let us now discuss the parameter Cq. For the systems with transla- 
tional invariance || |8| the parameter Cq was determined only in the special 
limit case, the so called hydrodynamic regime. As we see below, in both 
hydrodynamic and collisionless regimes the dissipation can be neglected and 
thus the Hamiltonian approach is valid. As a result the Eqs. (2.9-10) are 
valid in both regimes, but with different values of the parameter Cq. For 
example, as follows from Refs. @, §, in the systems with the translational 
invariance the collisionless value of the parameter Cq in Eqs. (2.9-10) is zero. 
But for superconductors with several electronic bands the situation is more 
complicated and the parameter Cq can be nonzero even in the collisionless 
regime. Moreover one can have the collisionless regime for one band and the 
hydrodynamic regime for the other band. 

3. Particles and holes contributions. The low-frequency dynamics and 
thermodynamics of the fermi-liquid or superconductors are determined by 
the low-energy quasiparticles. In the same way the low-frequency dynamics 
of vortices is determined by the low-energy excitations in the vortex core. 




k x [p(v s - v L ) + C (v L - v n ) + p n {y n - y s 



)] = 



(2.10) 
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The latter are concentrated on the anomalous branch of the spectrum [[11] . 
Here we follow the simplified version of [l], § (see Refs.||, |[]). Let us start 
with the axisymmetric vortex in the translational invariant surrounding. The 
spectrum of the low-energy excitations in the core is defined in the frame of 
the moving vortex, where the Hamiltonian does not depend on time, if there 
are impurities. The energy of the low-energy branch is expressed in terms 
of the canonically conjugated variables: the angular momentum Q and the 
angle a of the linear momentum, k = (hp cos 6, /c^ sin ^ cos a, /cFsin^sina), 
in the transverse plane: 

H = CuoQ + v s • k . (3.1) 
The first term describes the orbital motion of fermions bound in the core 



around the vortex axis ||11JI - The effect of the vortex on the motion of the 
quasiparticle is similar to the magnetic field. The quantization of this orbital 
motion leads to the descrete levels (see also Sec. 4) with either integer or 
half-odd integer generalized angular momentum Q ||12|| . The frequency of 
rotation u is the function of 9. The direction of rotation is determined by 
the "chirality" factor C = ±1. For the conventional case of the particle-like 
excitations in the vortex with winding number N — 1 one has C — — 1 [[11 



while C = 1 for the hole-like excitations in the vortex with the same winding 
number. This follows from the index theorem, which relates the number of 
fermion zero modes in the vortex core to the vortex winding number N [0]: 
in the case of holes the topological invariant, which determines the number 
of zero modes, changes sign. 

The second term is the energy due to the superflow in the vortex frame. 

When the vortex moves with respect to the heat bath, its dynamics is 
nonequilibrium. The kinetics of the fermions on the branch in Eq.(3.1) is 
governed by the Boltzmann equation for the distribution function, n(Q,a), 
in the Q — a phase space ||]: 

dn dn <9((v s - v L ) • k) dn n{Q , a) - n eq (Q , a) 

m +CuJ °d^ d~a dQ- r (3 ' 2) 

The last term describes the relaxation to the equilibrium distribution function 
n eq determined by the heat bath outside the core: 

n eq (Q, a) = f(H - v n ■ k) = f(Cuj Q + (v s - v n ) ■ k) , (3.3) 
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where f(E) = (1 + exp(E /T)) 1 is the Fermi-function. 

If r does not depend on Q one gets the equation for average momentum 

9 t k+C^ zxk+^4sin^zx(v n -v L )(/(A(T))-/(-A(T))) = -- , (3.4) 

4 T 



r doi 

dQ—(n(l,a)-n eq (Q,a))k . (3.5) 
J In 



.- 1 r ,„da 
k = - 
2 

Here we take into account that / dQdgn is limited by the bound states 
below the gap A(T), since above the gap A(T) the spectrum of fermions is 
continuous Q. The effective interlevel distance for the unbound (delocalized) 
states is ujq = and they will be considered below. 

In the steady state of the vortex motion one has <9tk = and the Eq.(3.4) 
is easily solved. The rhs of Eq.(3.4) gives the momentum flow to the heat 
bath and thus the following force due to bound states below A(T) 



loc 



I ^£^,2 j n 2 q tanh ^— -f(Vr — V n )LUnT — C& X (Vr — V n )l . 



dk Zl „ 2 „ ; , 2fl 1 ^ A(T), 

(3.6) 

The spectral flow of unbound states above A(T) is not suppressed, since the 
corresponding u t = 0. This gives 



deloc 



-C|£ y dcostf sin 2 # (l -tanh^p) z x (v L - v„) . (3.7) 



Thus the total nondissipative spectral-flow force is 
k 3 



nondiss sp flow 



-C-^- ( dcos6 sin 2 6 
8tt J 



lu 2 t 2 AT 

1 tanh -A-Z 

1 + uj 2 t 2 2T 



zx(v L -v n ) 
(3.8) 



The dissipative part of the spectral-flow force 

Fdiss sp flow = (v n - v L )-^ / dcos6» sin 2 9 tanh — — g-— , (3.9) 



S7T 



2T 1 + c^r 2 



can be neglected when A = u r/(l + uj 2 t 2 ) <C 1. So, the vortex motion is 
governed by the conservative Hamiltonian dynamics either in the so called 
hydrodynamic regime, when ujqt <C 1, or in collisional regime, when ujqt 3> 1. 



6 



In both cases the contribution of the spectral flow to the parameter Co in 
Eq.(2.10) from one spin direction acquires universal values at low T: 



C 

— = , lu t > 1 , (3.10) 
m 

C 

— = V P , ujqt <C 1 , particle states , (3- 11) 

m 

C 

— = -V h , cu t < 1 , hole states . (3.12) 
m 

Here V p = p 3 F /6n 2 is the volume of the Fermi-sphere of particle states (with 
one spin direction), while V h = pp/6ir 2 is the volume of the Fermi-sphere of 
the hole states. 

In conventional superconductors there are no zeroes of energy in the quasi- 
particle spectrum and thus the Fermi surface is absent. Zeroes however ap- 
pear due to vortices. In the continuous description of the vortex core, the 



zeroes in the classical energy spectrum E = ye 2 (k) + |A(k, r)| 2 are concen- 
trated on the vortex axis, but can be splitted into point nodes distributed 



in the vortex core [13]. The function k(r), which shows the position of zero 
of E in the momentum space as a function of the coordinate r in the real 
space, maps the cross section of the vortex into the momentum space. Thus, 
if one sweeps the cross section of the iV-quantum vortex, one obtains the 
closed surface in the momentum space, swept by zeroes. The volume within 
this surface is just Cq/ui times iV ||. The physical meaning of Cq/tti is the 
number of the electronic states which remain kept by the heat bath during 
the vortex motion. In the frame moving with the vortex this corresponds 
to the number of the states flowing from the vortex to the heat bath in the 
extreme spectral flow regime. 

4- Electronic bands in crystals, open orbits. Let us consider the crystal 
with anisotropic Fermi-surface. In the case of the particle states the Eq.(3.11) 
remains to be valid for the arbitrary closed surface of zeroes with V p being 
the volume within the surface ||. If one neglects the dependence of ujqt on 
a, 6 and Q one can write an interpolating equation for the closed surface of 
particle states: 

C (particles) V p 



m 1 + UJnT 2 



(4.1) 
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The contribution of the holes to Cq is 



= y B _ V *_L_ = yB^lL- + v, 1^ , (4.2) 
m 1 + cjgr 2 1 + u;qT 2 1 + UqT 1 

where V B is the total volume of the Brillouin zone, V h = V B — V p is the 
volume of the hole states. The Eq.(4.2) follows from two arguments. (1) 
The spectral flow parameter C changes sign for holes, as was discussed in 
the previous Section. (2) The completely filled band should be considered 
as uneffected by the motion of the vortex, the electrons on these bands are 
completely pinned by the crystal lattice. This corresponds to the limit of the 
extremely fast relaxation and thus to the extreme spectral flow. 

The Eq.(4.2) agrees with the microscopic calculations in Ref. 0. If the 
gap A is small compared to the Fermi energy, then p is very close to the 
volume of all the particle states: p ~ mJ2 V p . If there is only one band and 
this band contains holes, then p « m(V B — V h ). In the superclean regime, 
ojqt > 1, the Eq.(4.2) gives p - C ~ m(V B - V h ) - mV B = -mV h . As a 
result <7Haii = — (ec/B)V h in agreement with Ref.[Q]. 

The Eqs.(4.1) and (4.2) should transform into each other during the fill- 
ing of the Brillouin zone, when the Fermi surface of particles transforms to 
the Fermi surface of holes. However Eqs.(4.1) and (4.2) differ by the value 
V b uJqT 2 / (1 + uJqT 2 ). The key to the problem of matching these equations in 
the Lifshitz transition is provided by open orbits, which appear as intermedi- 
ate stage between particle and hole Fermi-surfaces. The open Fermi-surfaces 
with complicated topology were discussed in the relation to the Hall effect 



in the normal metal (see recent paper JTj] and references there). Here we 
show that in the presence of the open surfaces of zeroes the parameter ujqt 
is small and Eqs.(4.1) and (4.2) match each other. 

In the semiclassical approach the energy of fermions in the vortex core: 

£ 2 = e 2 (k) + |A(r)| 2 . (4.3) 

The lowest energy levels are concentrated in the vicinity of zeroes k of the 
spectrum e(k), ie close to the former Fermi-surface of the normal metal. Near 
the Fermi-surface one has e(k) = v(k ) ■ (k — k ) = —iv±d s , where s is the 
coordinate along vj_ = v(k ) — z(z ■ v(k )). Close to the vortex axis one has 

|A(r)| 2 ^ 7 2 r 2 = 7 2 (s 2 + 6 2 ) , (4.4) 
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where b is the impact parameter. We assume for simplicity that the core 
radius is large compared to the size of the electronic orbits: since the result 
is of the topological origin, it should not depend on the model. First we 
quantize the fast motion along s. According to supersymmetry the lowest 
energy level of this motion lies exactly at zero energy ||. As a result the 
spectrum of the excitations in the core is determined by the slow motion along 
b, ie along the line of the intersection of the Fermi-surface e(k) = with the 
plane k z = const. It is given by E(k\\,b) = 76, where k\\ is the coordinate 
along the line of zeroes in the momentum space. The quantization of the 
slow motion, § dk\\b(k») = 2im, with b(k») = E/'y gives the levels of bound 
states in the core 

1 1 rdku , . 

E n = -nuj , — = - f —± . (4.5) 

UJ 7 J Z7T 

For the closed spherical Fermi-surface this leads to the conventional result 
for the states in the vortex with large core radius: 

E n = -nu , u = ~, r- ^ • (4.6) 

kp sin 6 

For the open orbits the integral in Eq.(4.5) diverges which gives ujq = as it 
was expected. 

5. Discussion. The Eqs.(4.1) and (4.2) can be generalized to the case of 
several bands. If there are no open surfaces of zeroes, the total C^/m contains 
the positive contribution from the particles, the negative contribution from 
holes and the positive contribution kV B where k is the number of Brillouin 
zones which are either completely filled or contain the hole states: 

— = E Va-i — - E v b \ — + kyB ■ (s- 1 ) 

This also includes the summation over the spin indices. 

Each Fermi surface has its own ujqt. The conservative Hamiltonian ap- 
proach applies and the Eq.(5.1) holds, only if for each band the parameter 
A = co> t/(1 + WqT 2 ) < 1. If A are not small, the Eq.(5.1) can be considered 
only as interpolation since the real uq and r are complicated functions of the 
impact parameter and momentum p z . The limiting cases, when the Eq.(5.1) 
holds, can include the cases when uj q t is small in one zone and large in the 
other. Thus Cq as a function of external parameters (doping or direction of 
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magnetic field) should have plateaus interrupted by regions where one of the 
parameters A changes between and 1. The latter occurs also during the 
change of the topology of orbits. 

The Eq.(2.9) can be also applied to other inhomegeneous systems which 
become homogeneous on a large scale, such as Josephson junction arrays 
(J J A). In some cases the corresponding quantities p, p s and C can be ob- 
tained after averaging over the scale of the inhomogeneity. In the system of 
the SNS contacts, the parameter p — Co is again small in the hydrodynamic 



limit due to approximate particle- hole symmetry [T5|]. This leads to the al- 
most ballistic motion of vortices in the absence of the supercurrent, when 
j s = p s v s = 0. It is still unclear whether the approximate cancellation of p 
and Co occurs in SIS contacts. 
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